YITP-SB-02-25 

Toward Solving the Cosmological Constant Problem? 



Olindo Corradini^*, Alberto Iglesias^^ and Zurab Kakushadze^''' 
^ Dipartimento di Fisica, Universitd di Bologna and INFN, Sezione di Bologna 
via Irnerio 46, 1-40126 Bologna, Italy 
'^C.N. Yang Institute for Theoretical Physics 
State University of New York, Stony Brook, NY 11794 
^ Royal Bank of Canada Dominion Securities Corporation^ 
1 Liberty Plaza, 165 Broadway, New York, NY 10006 
(December 10, 2002) 



Abstract 



We discuss the cosmological constant problem in the context of higher codi- 
mension brane world scenarios with infinite- volume extra dimensions. In par- 
ticular, by adding higher curvature terms in the bulk action we are able to 
find smooth solutions with the property that the 4-dimensional part of the 
brane world-volume is flat for a range of positive values of the brane tension. 
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I. INTRODUCTION AND SUMMARY 



As was proposed in [1], one can reproduce four-dimensional gravity on a 3-brane in 6 
or higher dimensional infinite-volume bulk if one includes the Einstein-Hilbert term on the 
brane (the latter is induced at the quantum level if the brane matter is non-conformal) . Grav- 
ity then is almost completely localized on the brane with only ultra-light modes penetrating 
into the bulk, so that gravity is four-dimensional at distance scales up to an ultra-large 
cross-over scale Tc (beyond which gravity becomes higher dimensional), which can be larger 
than the present Hubble size. In particular, this is the case for codimension-2 and higher 
tensionless branes [1] as well as for codimension-2 non-zero tension branes [2] . This is also 
expected to be the case for codimension-3 and higher non-zero tension branes [3,4]. 

A careful analysis of linearized gravity in such backgrounds requires smoothing out higher 
codimension singularities [1,3]. This is already the case for tensionless branes, where the 
background is non-singular (in fact, it is fiat^), but the graviton propagator is singular. How- 
ever, as was pointed out in [4], these singularities in the graviton propagator are cured once 
we include higher curvature terms on the brane. These terms are induced at the quantum 
level along with the Einstcin-Hilbcrt term, and provide an ultra-violet cut-off for the gravi- 
ton propagator. In the presence of these terms, to reproduce four- dimensional gravity up 
to the present Hubble size, we need to assume that the bulk Planck scale Mp ^ (1 mm)~^ 
regardless of the codimension of the brane (assuming that it is 2 or higher) [4] . 

In the case of non-zero tension branes, however, the situation is more complicated as in 
the case of (5-function-like branes the background itself becomes singular (for phenomeno- 
logically interesting non-BPS branes on which we focus in this paper). More precisely, 
here we focus on solutions with vanishing cosmological constant on the brane. Then in the 
codimension-2 case the singularity is very mild as in the extra 2 dimensions the background 
has the form of a wedge with a deficit angle, so the singularity is a simple conical one [2]. 
As was discussed in [3], this singularity can be consistently smoothed out. In this case 
gravity on the brane was analyzed in [2,3], where it was found that the behavior of gravity 
is essentially unmodified compared with the tensionless brane cases. In particular, gravity 
is almost completely localized with only ultra-light modes penetrating into the bulk, and 
four-dimensional gravity on the brane can be reproduced up to distance scales of order of 
the current Hubble size. As to the brane tension, in codimension-2 solutions with vanishing 
brane cosmological constant it can take continuous values in the range bounded below by 
zero and bounded above by its critical value [2], which is of order Mp < (1 mm)~^ (for a 
3-brane) . Since this does not improve the current experimental bound on the 4-dimensional 
cosmological constant, we are prompted to consider codimcnsion-3 and higher solutions. 

In codimension-3 and higher cases with non-zero tension branes the singularities in the 
background itself are more severe. Actually, there are two types of singularities we must 
distinguish here. Thus, we can have a singularity at the origin r = of the extra space, 
that is, at the location of the brane. These singularities can sometimes be simple coordinate 
singularities, in which case they are harmless. On the other hand, true r — singularities 



More precisely, in the absence of higher curvature terms in the bulk the latter can be Ricci-flat 
(here and in the following we focus on the solutions with vanishing brane cosmological constant). 



2 



can be smoothed out using a procedure discussed in [3]. This procedure goes as fohows. 
Consider a codimension-d {d > 2) 5-function-hke source brane in /^-dimensional bulk. Let 
the bulk action simply be the D-dimensional Einstein-Hilbert action, while on the brane 
we have the induced {D — d) -dimensional Einstein-Hilbert term as well as the cosmological 
term corresponding to the brane tension. As we have already mentioned, the background 
in this case is singular [5]. One way to smooth out such a singularity is to replace the 
{D — (i)-dimensional world-volume of the brane by its product with a rf-dimensional ball 
of some non-zero radius e. As was pointed out in [3], in this case already for a tensionless 
brane the gravitational modes on the brane contain an infinite tower of tachyonic modes. 
This can be circumvented by considering a partial smoothing out where one replaces the 
(D — d) -dimensional world-volume of the brane by its product with a {d— l)-sphere S"^"^ of 
radius e [3]. As was pointed out in [3], this suffices for smoothing out higher codimension 
singularities in the graviton propagator as in the codimension- 1 case the propagator is non- 
singular [6]. Moreover, in the case of tensionless branes as well as in the case of a non-zero 
tension codimension-2 brane we then have only one tachyonic mode (with ultra- low — m^) 
which is expected to be an artifact of not including non-local operators on the brane. As 
to the background itself, this procedure does smooth out r = singularities. However, not 
surprisingly, it does not suffice for smoothing out the second type of singularities [7] (such 
singularities are absent in the codimension-2 case but are present in higher codimension 
cases with non-zero brane tension), which we discuss next. 

Thus, in codimension-3 and higher cases with vanishing brane cosmological constant we 
can have singularities at r = r* with some non- vanishing r*, which depends on the brane 
tension [5]. The aforementioned smoothing out procedure does not cure such singularities. 
In [7] it was suggested that higher curvature terms in the bulk might be responsible for curing 
such singularities. The purpose of this paper is precisely to study effects of higher curvature 
terms in the bulk on such singularities. 

Not surprisingly, studying such backgrounds in the presence of higher curvature terms in 
the bulk is rather non-trivial. To make the problem more tractable, we focus on a special kind 
of higher curvature terms in the bulk. In particular, we consider adding the second-order (in 
curvature) Gauss-Bonnet combination in the bulk. This Gauss-Bonnet combination is an 
Euler invariant in four dimensions. In higher dimensions, however, it has a non-trivial effect 
on the equations of motion if the background is not flat. A simplifying feature of the Gauss- 
Bonnet combination is that it does not introduce terms with third and fourth derivatives in 
the background equations of motions, but rather makes them even more non-linear. Albeit 
non-trivial, these equations can in certain cases be analyzed analytically, so we can get some 
insight into the effect of higher curvature terms on the aforementioned r — singularities. 

The results of our analyses, which we present in the remainder of this paper, suggest that 
higher curvature terms in the bulk might indeed cure these singularities. To simplify our 
discussion, as far as explicit computations are concerned, we focus on codimension-3 cases, 
where we have additional simplifications in the equations of motion (higher codimension 
cases are not expected to be qualitatively different). In particular, we explicitly study 
examples of: a string in 5D bulk, a membrane in 6D bulk as well as a 3-brane in 7D bulk. 
In these cases we argue that smooth solutions with vanishing brane cosmological constant 
do exist for a continuous range of positive brane tension. In the remainder of this paper, 
following the terminology adopted in [8] , we refer to such solutions as diluting solutions. 
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Before we turn to a detailed description of our solutions, let us describe the geometry 
of the aforementioned diluting solution corresponding to a 3-brane in 7D bulk. In the extra 
three dimensions we have a radially symmetric solution where a 2-sphere is fibered over a 
semi-infinite line [ro, oo). The space is curved for tq < r < e, at r = e we have a jump in the 
radial derivatives of the warp factors (because r = e is where the brane tension is localized) , 
and for r > e the space is actually flat. So an outside observer located at r > e thinks that 
the brane is tensionless, while an observer inside of the sphere, that is, at r < e sees the 
space highly curved. It would take this observer infinite time to reach the point r = Tq, 
where we have a coordinate singularity (that is, the corresponding geodesies are complete, 
so this is not a naked singularity but a coordinate one). This is an important point, in 
particular, we did not find smooth solutions where the space would be curved outside but 
fiat inside of the sphere. And this is just as well. Indeed, as was shown in [7], if we have 
only the Einstein-Hilbert term in the bulk, then we have no smooth solutions whatsoever 
(that is, smoothing out the 3-brane by making it into a 5-brane with two dimensions curled 
up into a 2-sphere does not suffice) . What is different in our solutions is that we have added 
higher curvature terms in the bulk, which we expected to smooth out some singularities. 
But serving as an ultra-violet cut-off higher curvature terms could only possibly smooth out 
a real singularity at r < e, not at r > e. And this is precisely what happens in our solution 
- the presence of higher curvature terms ensures that we have only a coordinate singularity 
at r < e instead of a real naked singularity as would be the case had we included only the 
Einstein-Hilbert term. 



II. SETUP 

The brane world model we study in this paper is described by the following action: 



R-A 



+ 



R + ^ ( ^RmN + ^MNPR 



(1) 



Here Mp is the (reduced) D— dimensional Planck mass, while Mp is the (reduced) {D — 1)- 
dimensional Planck mass; E is a source brane, whose geometry is given by the product 
YiD-d-1,1 y. s^~i, where R^-'^-i'i is the {D — d) -dimensional Minkowski space, and is a 
{d— l)-sphere of radius e (in the following we will assume that 0? > 3). The quantity Mp'^A 
plays the role of the tension of the brane E. Also, 



MN 



(2) 



where x"^ are the {D — 1) coordinates along the brane (the D-dimensional coordinates are 
given by x^ — {x"^,r), where r > is a non-compact radial coordinate transverse to the 
brane, and the signature of the D-dimensional metric is (— , +,..., +)); finally, the {D — 1)- 

dimensional Ricci scalar R is constructed from the {D — l)-dimensional metric G^^y. In the 
following we will use the notation = (a;",r), where x°' are the {d— 1) angular coordinates 
on the sphere Sf~^. Moreover, the metric for the coordinates will be (conformally) fiat: 



5ij dx^dx^ 



dr^ _l_ r^'j^p dx^dx^ 



(3) 
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where 7„a is the metric on a unit (rf— l)-sphere. Also, we will denote the {D — d) Minkowski 
coordinates on R,'^^'^^^'^ via (note that x"* = {x'^,x°')). 
The equations of motion read 



Rmn — -^GmnR + i 



+ 



-G 



where S = Sm'^, and 



om f at 



Rmn r\Gmn i^R 



L 5(r - e) = , 



(4) 



— I? r? o t? t?'S' _l t? t? i^ST n tdRS d 

= ^-tlMN — ^-n-MS-ti- N + -n-MRST-n-N — ^-K -n-MRNS 

D-3 /]\/rD-2 



(5) 
(6) 



Here we are interested in solutions with vanishing {D — d) -dimensional cosmological con- 
stant, which, at the same time, are radially symmetric in the extra d dimensions. The 
corresponding ansatz for the background metric reads: 



ds'^ = exp(2A) rjn,, dx^dx" -\- exp(2S) 5ij dx^'dx^ , 



(7) 



where A and B are functions of r but are independent of x'^ and x". We then have (here 
prime denotes derivative w.r.t. r): 



R/xi^ — , 
R^{d-1) X, 

R^. = -Vi^r^e^^^'^^ U" + id- 1)-A' +iD- d)iA'f + id- 2)A'B' 



Rj.j. — — (c? — 1) 
R = 



+ {D-d) \a'B' - {Af - A" 



B" + -B' 
r 

B" + {2d - 3)^5' + {d- 2){Bf + {D- d)A'B' + (D - d)^A' 

2(d - 1)B" + 2(d - 1?-B' + 2(D - d)A" + 2(D - d)id - 1)-A' + 

r r 



{d -l){d- 2) {B'Y + {D-d){D-d+l) (A'y + 2{D-d){d- 2)A'B 



where 



A = 



d-2 



-2B(e) 



(8) 
(9) 

(10) 
(11) 
(12) 

(13) 
(14) 

(15) 



The equations of motion then read: 

)-d{d - l){Af + d{d - 1)N + hd -l){d- 2)S ■ 



fhfl 



e e-'^|-d(d -l){d- 2){d - 3){A'y + 2d{d -l){d- l){d - 2)N' + 
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1 



-{d -l){d- 2){d -3){d- 4)S' + 2d{d -l){d- 2){d - l){A'yN 



+d{d -l){d- l){d - 2){A'yS + 2d{d -l){d- 2){d - 2>)NS | = , 
]-(d -l)(d- 2){A'f + (d- 1)M +(d-l){d- 1)N +{d- 1)P + hd -l){d- 2)S - 
i e-'^|^(d -l){d- 2){d - 3){d - A){A'f + 2{d - l){d - 2){d -l){d- 2)N^ + 

^{d -l){d- 2){d -3){d- 4)^2 + 2(d -l)(d- 2){d - ?>){A'fM + 

2{d- l)(d-2){d-?>){d-l){A'fN + 2{d-l)(d-2){d~l){A'fP + 
(d. -l)id- 2){d - - 2){AyS + A{d -l)(d- 2){d - l)MN + 
2(d- l){d-l){d-2)MS + A{d-l){d-l){d-2)NP + 

2{d-l){d-l){d-2){d-2>)NS + 2{d-l){d-2){d-2,)PS^ 



(16) 



+ 



A - (d-l)A 



\d(d - l){A'f + dM + d{d - 2)N +{d-2)P+hd- 2){d - 3)S - 
i e-2^(^d(d -\){d- 2){d - 3){A'y + 2d{d -l){d- 2){d - 3)N^ + 



(17) 



1 



l2 

{d -2){d- 3){d -4){d- 5)S^ + 2d{d -l){d- 2){A'fM + 



2d{d -l){d- 2){d - 2){A'YN + 2d{d - l){d - 2)(A')'P + 

d{d -l){d- 2){d - 3){A'fS + 4d(d -l){d- 2)MN + 2d{d -2){d- 3)MS + 

4d{d -2){d- 3)NP + 2d{d -2){d- 3){d - A)NS + 2{d -2){d- 3){d - A)PS^ 

1 



+ 



A - (rf-3)A 



L5{r - e) = , 



(18) 



where we have defined 



M = A' + {Af - A'B' , 
B' 



N = A'{B' + - 
r 



P = B" + 

r 

d = D-d . 



S = {B'Y + 2 



B' 



Above the third equation is the (a/?) equation, the second equation is the (yuz/) equation, 
while the first equation is the (rr) equation. Note that the latter does not contain second 
derivatives of A and B. Also note that the other two equations do not contain third and 
fourth derivatives of A and B (this is a special property of the Gauss-Bonnet combination 
we mentioned in Introduction). 

The equations of motion (16)-(18) are highly non-linear and difficult to solve in the 
general case. However, in the d = A and especially d = 3 cases various (but not all) terms 
proportional to ^ vanish. This is due to the fact that the Gauss-Bonnet combination is an 
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Euler invariant in four dimensions. To make our task more tractable, from now on wc will 
focus on the codimension-3 case {d — 3). We do not expect higher codimension cases to be 
qualitatively different. 

Not only is the complexity of the above equations of motion sensitive to the value of d, 
but also to the value of d. In particular, we have substantial simplifications in the cases 
of d = 2 and d = 3 corresponding to the non-compact part of the brane S being a string 
respectively a membrane. Note that these simplifications are specific to the Gauss-Bonnet 
combination. In particular, if we set the Gauss-Bonnet coupling ^ to zero (that is, if we 
keep only the Einstein- Hilbert term in the bulk), there is nothing special about the d — 2,3 
cases. This suggest that our conclusions derived from explicit analytical computations for 
the d = 2,3 cases can be expected to hold in o? > 4 cases as well (in particular, in the case 
of a 3-brane in 7D and, as we mentioned above, even higher dimensional bulk). In fact, our 
analytical and numerical results in the case of a 3-brane in 7D bulk confirm this expectation. 

III. NO EINSTEIN-HILBERT TERM IN THE BULK 

To warm up and get a feeling for whether we can expect to find smooth solutions in the 
case where we have the Einstein-Hilbert as well as the Gauss-Bonnet terms in the bulk, in 
this section we will discuss a somewhat simpler (and more of a toy) problem. In particular, 
we will study the case where wc have no Einstein-Hilbert term in the bulk but only the 
Gauss-Bonnet term. This simplifies the equations of motion enough so that we can solve 
them analytically. In the next section we will come back to our original problem where we 
have both the Einstein-Hilbert term as well as the Gauss-Bonnet term in the bulk. 

A. String in 5D Bulk 

Let us first consider a string {d = 2) propagating in 5D bulk {d = 3). Our equations of 
motion then read: 



/\2 



(A') 



B' 2 



3(Sy + 6— + 



= 0, 



A" + {A'f - A'B' 
1 



(B'f + 2— 
r 



e''^{K-2)\L5{r-€) , 



(19) 



+ 2A' [B' + -\{b" + — 
rj\ r , 




(20) 



(21) 



where we have introduced the notation 

K = K/\. 

Let us now solve the above equations of motion. 



(22) 
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Thus, we have two types of solutions. 
• Solution A: 



r < e : A{r) — constant, B{r) — anything , 
r > e : A{r) = In (^^)" ' , B{r) = In (^yj 



Solution B: 



r<e: ^(r) = In f — V \ B{r) ^ln(—) 



\ r J 



\ r J 



Here 



r > e : A{r) — constant, B{r) — anything 



a — a± = 1 ± —z=z . 

V3 



(23) 
(24) 

(25) 
(26) 

(27) 



The jump conditions at r = e are given by: 

A' ({B'f + 2— ' " 



{A'y{B'+-^ 



_ K-2L 



.S(e) 



(28) 

(29) 



Note that, without losing any solutions parametrically, we can assume that B{r) = constant 
whenever A{r) — constant. The jump conditions then read: 



± 
± 



A' f(S')' + 2— 



{Ar{B' + -^ 



r=e± 



K-2L 



Me) 



r=e± 



8^6^" 



(30) 
(31) 



where the upper (+) sign corresponds to Solution A, while the lower (— ) sign corresponds 
to Solution B. We then have: 



± 



a{a — l){a — 2) 



e 8^ \ e J 



(32) 
(33) 



which gives the following solution for k: 

K^2(a- if . (34) 
The brane tension, in terms of the parameters of the solution, is therefore given by: 



A = 2 



-)■ 



(35) 
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Note that the brane tension is always positive (assuming that the "brane-width" e is non- 
zero). At first it might seem that it can be arbitrary as we can adjust the integration 
constant tb- However, we also have a relation between the Gauss-Bonnet coupling ^ and 



which implies that 



A=^. (37) 



That is, these solutions are not "diluting" , rather they exist only if we fine-tune the brane 
tension and the Gauss-Bonnet coupling (we discuss the reason for this in the next subsec- 
tion). In the following we will see that this specific to the case at hand, and diluting solutions 
do exist in other cases. 

Before we end this subsection, let us analyze the singularity structure in the above 
solutions. Singularities can potentially occur at r — > oo in Solution A and at r — > in 
Solution B. Thus, the line element in the non-trivially warped part of the space-time is 
given by: 

ds^ = (^^^ rj^^dx^dx" + (^y ^ Sijdx'dx^ , (38) 

which is singular at r = for both roots a = a±. However, only for a = a_ is the space 
truly singular, whereas for a — a+ we merely have a coordinate singularity. Thus, consider 
a 2n-derivative scalar. Such an object - let us call it x^"^ - has the following expression in 
terms of the extra-space radius 

^(n) ^ - ^2n(a-l)_ ^gg^ 

The latter blows up for a = q;_ as r approaches zero; this singularity is a naked singularity. 
One can indeed consider radial null geodesies with affine parameter a and use that Goodt/da 
is constant along geodesies to obtain 

^ ~ r^^-i . (40) 
da 

For a — these geodesies terminate with finite affine parameter as r approaches zero, 

a ~ r^(^"") + constant , (41) 

Thus, we have incomplete geodesies reaching a point of divergent curvature. On the other 
hand, for a = the expression (39) vanishes as r approaches zero, the aforementioned 
geodesies are complete {i.e., a — >• oo), and it is not difficult to see, by doing a similar calcu- 
lation, that radial time-like geodesies extend to infinite proper time in this limit. Therefore, 
r = is a coordinate singularity in this case. This then implies that for Solution B we must 
choose a — It then follows that the Gauss-Bonnet coupling ^ is positive in this case. 
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Similar considerations apply to the r — > oo singularity. In this case we have a naked 
singularity for a = while for a = we merely have a coordinate singularity. This 
implies that for Solution A we must choose a = Note that in this case the Gauss- 
Bonnet coupling ^ is also positive. 

Another comment we would like to make here concerns the volume of the extra space, 
which is infinite. Indeed, this volume is given by 

oo — 

dr r'^-^ e'^^ e'^^ . (42) 



In Solution B this volume is clearly infinite. In fact, in Solution B even the volume of the 
space inside of the {d — l)-sphere (that is, for < r < e) is infinite: 

On the other hand, in Solution A the volume of the space inside of the {d — l)-sphere is 
finite; but the volume outside of the {d — 1) -sphere (that is, for r > e) is infinite: 

(^^)*""" (^)"- ~ r-'-* . (44) 

That is, in both Solution A and Solution B the extra space has infinite volume^. 

Thus, as we see, we have sensible infinite-volume non-singular solutions if we take the 
bulk action to be given by the Gauss-Bonnet combination. However, as we have already 
mentioned, these solutions exist only for the fine-tuned value of the brane tension. If the 
brane tension is not fine-tuned, then we expect to have solutions where the non-compact 
part of the brane is infiating. 



B. Membrane in 6D Bulk 

Let us now consider the case of a membrane {d = 3) in 6D bulk {d = 3). As in the string 
case discussed in the previous subsection a drastic simplification of the equations of motion 
occurs if we neglect the Einstein-Hilbert term in the bulk action. In this subsection we will 
focus on this case. The (rr), (rr — a/3) and (iiu) equations of motion read: 

A'(2A'^,s{lJ-l,yo, (45) 

2N {N-M) + {A'f (S-P + N-M)^ -^e^^AL5{r - e) , (46) 

^4^ 

2N{N + 2P) + {A'f{2P + S) + 2M{2N + S) = ^e^^ (A - 2A) LS{r - e) . (47) 



^Notc that, in Solution A, the extra space would appear to have finite volume had we chosen 
a = a+. However, as we saw above, in this case we would have a naked singularity at r — > oo, 
and the finiteness of the volume of the extra space would be due to truncating the space, which 
is geodesically incomplete, at the naked singularity. A similar situation will arise again in the 
following sections. 
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Here 

p = r . 

If we use p instead of r, then the metric takes the following form: 

ds^ = exp(2A) 1]^^ dx^dx" + exp(2C) dp'^ + p^ 7^/3 dx'^dx'^ , 



(48) 



(49) 



where, as before, 7q,^ is the metric on a unit {d — l)-sphere, while exp(C) = exp{B) / p' . 

For our purposes here it will be convenient to use U = exp(— C) = r p'/p. Then from 
(45) and (46) we have (for convenience here we also include the equation for B): 



A' 



1 - 3C/^) , 



2rU 



r{9U^ + 2[/2 + 1)U' + -{-33U^ + 5U^ + - 1) 



^ = __L,3e3-AZ5(,_,). 



(50) 
(51) 
(52) 



The last equation can be integrated for r 7^ e. Thus, we have the following solution: 



In (—] = — 2-\/3 arctanh (V3U{r)^ + — arctan 



where Tq is an integration constant, and 



Vjr) 



1 , fU(r) 
— arctanh 

V- \ U- 



(53) 



u± = 



v± 



'±1 + 2^3 
11 

±11 + 22^3 
T8 + 6v^ 



Note that u^u_ = 1/vTT, and v+V- = vTT/4. 

Next, note that in both the r ^ as well as the r 
dominated by the last term, so in these limits we have 



In 



arctanh 



U- 



(54) 
(55) 

00 limits the r.h.s of (53) is 

(56) 



and 



U (r) w i-u- 



1-2 1 7)'" 



(57) 



where the "+" sign corresponds to the r ^ 00 limit, whereas the "— " sign corresponds 
to the r — > limit. In particular, the 2n-derivative scalar x^'^^ introduced in the previous 
subsection reads 



0, 



(58) 
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and, therefore, the solution imphcitly given in (53) is smooth at r = and r — > oo (since 
the corresponding singularities are coordinate singularities). In the above sign conventions 
we have the following asymptotic behavior for the warp factors. 



1-3x^2 



B ~ {±u^ - l)ln (—\ . (60) 

Thus, we can check that null and time-like radial geodesies are complete in these limits. The 
affine parameter for null geodesies diverges as 



± 



(7 ~ r V / — > oo , (61) 

while the proper time diverges as r ~ r^"" — > oo. 

What about finite distance singularities? First, note that for the function U{r) in (53) 
we have 

^ 11 {?>u^-i){u'-vr){u^ + ui) 

Since 3rt^ < 1, the l.h.s. of this equation is strictly positive for U between — it_ and 
This together with (57) implies that U{r) is a monotonically increasing function bounded 
from below by its value —u_ at r = and bounded from above by its asymptotic value +u_ 
at r — s> oo. In particular, U{r) is finite. However, U{r) does go through zero at a finite value 
of r. According to (50) and (51) this could potentially lead to to a nasty finite-distance 
singularity in the full solution^. However, as we will now see, this does not actually take 
place for a range of values of the brane tension. 

Thus, we need to accommodate the jump conditions 



1 . . „.,o.o.. „.,o. ^ I g3B(.)^^ ^ (g3) 



_(_l + 3[/2)2(l + 3^2 



1 



-(l + t/')(-l + 3[/2) 



-e=^^(^)(A-2A)Z, (64) 



which follow from (46) and (47). To do this, we can consider two types of solutions: 

• Solution A: 

r < e : A = constant , B = constant {U = 1) . (65) 

• Solution B: 

r > e : A — constant , B — constant {U — 1) . (66) 



^As r — ro, the curvature scalar diverges as i? ~ (r — ro) ^, and radial null geodesies are 
incomplete. 
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We must check that we can consistently glue these solutions with the non-trivial solution 
with U{r) given by (53) for r > e for Solution A and for r < e for Solution B. 
Let us start with Solution B. The jump condition (63) then reads: 



^(-l + 3C/2)2(l + 3[/2) 



16 - ^ e^(^)Z . (67) 



Now consider solutions with U{e—) < (but U{e—) > —uJ). From our previous discussion 
it follows that such a solution is non-singular. Assuming that A > 0, this requires that ^ 
be positive - indeed, the l.h.s. of (67) is negative. We can then choose B{e) such that (67) 
is satisfied (note that U does not depend on B{e) since it is only related to B'). Next, the 
second jump condition (64) reads: 



^(1 + [/')(-! + 3C/') 



4_i^ll-^e^(^)Z. (68) 
4^ 



The l.h.s. of this equation is positive. This condition then can be satisfied if we assume that 
n <2 (note, however, that this condition is not necessary"^). In fact, this solution is diluting, 
that is, there is a continuous range of positive values of A for which this solution exists. To 
see this, consider the case where |C/(e— )| <^ 1. It is then not difficult to see that we have 

K « 2 , (69) 



This implies that 



K = K\^U\e-) . (71) 



Note the difference between this result and the corresponding result (37) for the string 
discussed in the previous subsection. In the latter case the brane tension had to be fine- 
tuned. In the case of the membrane we actually have diluting - the solution exists for a 
range of the brane tension A. Indeed, U{e—) depends on the integration constant ro, which 
is the source of this diluting property. Thus, in the case where |6'^(e— )| <^ 1 we have 

U{e-) ^ -I In i'-^] , (72) 



2 V e 

so that the brane tension A depends on tq. Note that in the case of the string the function 
U (r) defined as above is actually a numerical constant away from the brane. This is precisely 
why we did not find a diluting solution in the case of the string. On the other hand, for a 
membrane as well as higher dimensional branes this function is non-trivial and depends on 



^In fact, K > 2 values are "phenomenologically" most interesting - see the end of this section. 
However, since here we are discussing a membrane (and not a 3-brane), we will be somewhat 
cavalier. 
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an additional integration constant, hence the diluting property of such solutions. In turn, 
the reason why this is possible, say, in the membrane case is that equations of motion are 
more non-linear than those in the string case. 

Let us now consider Solution A. In this case the jump conditions read: 



1 

IP 



(-1 + 3C/2)2(1 + 
l + t/2)(-l + 3[/2) 



4 + 



4e 



(73) 
(74) 



If we take U{t+) > (but U{e+) < U-), then we have a smooth solution. In fact, this 
solution is also diluting. This can be seen by noting that choosing positive ^ and k < 2 
ensures that the jump conditions can be satisfied. 

The last comment we would like to make here concerns the volume of the extra space. 
It is evident that this volume is infinite in Solution B. For solution A we have the following 
expression for this volume: 



f 

Jo 



dr e^^ e^^ . 



Asymptotically at r — > oo we have 



B{r) ~ (ii_ - 1) In (—] . 

\rB/ 



(75) 

(76) 
(77) 



It is then not difficult to see that the volume of the extra space in Solution A is infinite 
as well. Numerical results (obtained using Maple 6) for A{r) and B{r) for Solution A are 
presented in Fig.l and Fig. 2, respectively. 

Thus, as we see, we have sensible infinite- volume non-singular solutions in the case of 
the membrane as well. Unlike in the string case, however, these solutions have the diluting 
property - they exist for a continuous range of the brane tension values. 



C. 3-Brane in 7D Bulk 

In this subsection we consider the case of a 3-brane {d — 4) in 7D bulk (d = 3). As in 
the previous subsections we will neglect the Einstein-Hilbert term in the bulk action. Let 



U = rB' + 1 . 
Then the (rr) equation of motion reads: 



(^0 



/\2 



{A'Y + 8^' - + 2 



, U ^ 3C/2 - 1 



which implies that either A' = or 

rA' = -AU ± VlOC/2 + 2 = -AU ± T 



0, 



(78) 



(79) 



(80) 



14 



The (rr — a(5) and {^v) equations of motion read: 
2[M-N] UAf + iVl + {Af [P - 5] = -L e'^AL(5(r - e) , 



(81) 



{Af + n\[P + N] + mUn + S + {Af] + {AfS = e^^ (A - 2A) LS{r - e) . (82) 



Note that 



r 

^- A'U^\ \r)TU - AU^ 



N 



U' 1 r 

M = — [10r]U - 4T] + — 30U^ - 9r]TU + 2 



(83) 
(84) 
(85) 
(86) 



where 77 = ±1. Using these expressions we can rewrite (81) as follows: 

12r;- 



{rjT - AU) \ rU' 



T 



AOOU'' - 127r]TU' + 56U - 5r]T 



(87) 



Thus, away from the brane the function U satisfies the following first order differential 
equation (assuming a non-trivial solution, that is ^4 7^ const.): 



rU' = 



mU'^ - 127r/Tf/2 + 56U - 5r]T 
40U - 13r]T+12r]/T 



(88) 



It is not difficult to show that for 77 = +1 the denominator is always negative, while for 
?7 = — 1 it is always positive. As to the numerator, it is positive for U < Ui and U > U2, it 
is negative for ui < U < U2, and it vanishes for U = Ui and U = U2, where 



T] — —1 : Ui — Ui ^ —.58 , U2 — U2 ^ —.43 , 
77 = +1 : ui — uf^ +.43 , U2—U2~ +.58 . 



(89) 
(90) 



This implies that we have the following solutions for U{r), all of which are monotonically 
increasing: 



7] = —1 
T] = —1 

r] = +1 



U{r ro+) —00 , U{r ^ 00) , 
f/(r — i> 0) — s> M2 , U{r- ^ ^Q—) +00 , 
U{r ^ 0) ^ < , U{r ^oo)^ut . 



(91) 
(92) 
(93) 



Here tq in each of the rj = —1 cases is an integration constant (see below). Let us discuss 
these solutions in more detail. 
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The 77 = +1 Solutions 



In this solution the r dependence of U is given by 

U{r) = /(ln(r/ro)) , 



(94) 



where tq is an integration constant, while the function f[x) is the solution of the differential 
equation 



d[_ _ Amp - mWp + 56/ - 5W 
dx ~ 40/-13iy + 12/H/ 



(95) 

00) = Ui and 



where W = a/10/^ + 2, and the boundary conditions are chosen as /(a: - 
f{x — > +00) = Note that f{x) is a monotonically increasing function. 

Next, to find a full solution, we must impose the jump conditions. These are given by: 



■46 


rp3 


290 


.15 


3 


■13 


rpS 


40 


.30 


— 
3 



24 

U^-UU- — T 



48^' ' 



6 1 ^+ 
U^--T 
5 



24^ 



(96) 
(97) 



Let us consider solutions where A and B are constant for r < e (in this region we then have 
U = 1), and have non-trivial profiles for r > e (in this region U is non-trivial). The jump 
conditions then read: 



■46 


rp3 


290 


.15 


3 


■13 






.30 




f' 



o 24 ■ 

U^-UU T 

5 



6 
5 



8^3- 



40 



2AC 



Let us rewrite these conditions as follows 
290 



24 ■ 
UU T 



46^3 

15 3 5 

57 o , 21 ■ 

— -90 -UU T 

20 5 



^64v^- — + — e^(^)Z 
.+ 3 48e 

= 60^3 - 104 + e^^'^L . 

e+ 24^ 



(98) 
(99) 

(100) 
(101) 



Then one can check numerically that we do not have a solution with positive k (that is, 
for positive A). Thus, k ranges between about —1.60 and —4.80 depending on the value of 
U{e+) (which must be between uf and Wg), while eexp[i?(e)] ranges between about 4.45 
and 1.85 times ^/Z, so that A ranges between about —.080 and —1.41 times Z^/^^. That is, 
this solution is diluting. Note that in this solution ^ must be positive. Also, the volume of 
the extra space in this solution is infinite. Thus, the volume of the extra space is given by: 



f 

Jo 



(102) 



On the other hand, for large r we have 
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A{t) ~ (^10(«+)2 + 2 - ^ut) In ( - ) , (103) 

B{r) ~ {ut - 1) In (^-^^ . (104) 

Then it is not difficult to show that the volume of the extra space is infinite, albeit the brane 
tension must be negative. 

Alternatively, we can consider solutions where A and B are constant for r > e. Note 
that in these solutions the volume of the extra space is automatically infinite. The matching 
conditions now read: 



46 „ 290 , „ 24 ^ 

— -UU T 

15 3 5 . 

57 o o 21 ■ 

— T^- 90 U^-UU T 

20 5 



^64V3- — e^(^)Z, (105) 
3 48^ ' ^ ' 



= 60^3 - 104 - e^^'^L . (106) 
24^ 



It is not difficult to show that in these solutions, which are diluting with negative brane 
tension, ^ must be negative. 

The T) — —1 Solutions 

Let us now discuss the rj — —1 solutions. If we choose the solution such that U{r ^ 
oo) — > Ui, then the r dependence of U is given by 

U{r) = /(ln(r/ro)) , (107) 

where tq is an integration constant, while the function f[x) is the solution of the differential 
equation 

4f ^ 400/^ + mWp + 56/ + 

dx 40/ + 13VF + 12/V^ ' ^ ' 

where W = y/lOp + 2, and the boundary conditions are chosen as f{x — > 0+) — > — oo and 
f(x — > +oo) = Ui- Note that f{x) is a monotonically increasing function. To ensure the 
absence of a singularity at r = rg (see below) we can consider the solution where A and B 
are constant for r < e, where e > Tq. However, it is not difficult to show that the volume of 
the extra space in this solution is finite. 

A solution with infinite- volume extra space can be obtained if we take A and B to be 
constant for r > e. In this solution the r dependence of U is given by 

U{r) = /(ln(r/ro)) , (109) 

where tq is an integration constant, while the function f{x) is the solution of the differential 
equation 

dl ^ 400/3 + mWp + 56/ + 5W 

dx 4Qf + 13W + 12/W ' ^ ^ 

where W = \/lOP + 2, and the boundary conditions are chosen as /(x — > — oo) = U2 
and f{x 0—) = +00. Note that f{x) is a monotonically increasing function. To ensure 
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the absence of a singularity at r = ro (see below) we must choose e < tq. The matching 
conditions then read: 



46 o 290 , o 24 ^ 

— - — - UU + — T 

15 3 5 . 

57 o o 21 

■— -90U^ -UU + —T 

20 5 



= _64V3-f -^e-(.)L. (Ill) 



= -60^3 - 104 - ^ e^^'^L . (112) 

In this solution must be negative for t/(e— ) < 1, while for U{e—) > 1 it must be positive 
(note that U{e—) in this solution can take values larger than Ug)- That is, we have two 
different branches here. The parameter k is not very sensitive to the value of C/(e— ), in 
particular, it smoothly changes from the ^ < branch to the ^ > branch - it ranges 
between about 2.13 if U{e—) is close to and 2.15 if U{e—) ^ 1. As to A, if U{e—) < 1, 
it grows from about 8.57 x 10~^ (L^/^^) if U{e—) is close to U2 up to infinity as U{e—) 
approaches 1. On the other hand, if U{e—) > 1, A decreases from infinity, as U{e—) moves 
away from 1, down to zero. Thus, we have a diluting solution with positive brane tension. 

Singularity Structure 

Before we end this subsection, let us analyze the singularity structure of the above 
solutions. In particular, ii U ^ U^, in any given solution for r — > or r — > 00, then the 
2n-derivative scalar introduced in the previous subsections behaves as 

^(n) ^ ~ , (113) 

while the afSne parameter a and proper time r, for null and time-like geodesies, respectively, 
have the following asymptotic expressions for the rj — ±1 solutions: 



^ _ ^-3U,±^iou^+2 ^ constant , (114) 
T ~ r^* + constant . (115) 

In the 77 = +1 solution U^, = for r — > and = U2 for r ^ 00. We therefore would 
have a coordinate singularity at r ^ 00 as x'-"^ — > and A and r diverge, but at r — > 
we would have a true singularity as we have diverging curvature and incomplete geodesies. 
This imphes that for 77 = +1 we must choose the solution where A and B are constant for 
r < e. Recall that in this solution ^ must be positive (as the solution with negative ^ has a 
true singularity at r = 0). 

In the rj = —1 case the solution with U{r 00) Ui has a true singularity at r — 00. 
This explains why in this solution we found that the volume of the extra space is finite. This 
appears to be an artifact of cutting off the extra space at the singularity. At any rate, since 
this solution is truly singular, we will no longer focus on it. 

On the other hand, the 77 = —1 solution with U{r — 0) — * U2 only has a coordinate 
singularity at r = 0, so this solution is perfectly acceptable. In this solution, which is diluting 
with positive brane tension, the volume of the extra space is infinite. Numerical results for 
this solution are given in Fig. 3 and Fig. 4. 

Finally, let us show that at r — > ro± in the the 77 = —1 solutions where C/(r — > rQ+) — > 
—00 and U{r ^ Tq—) — > +00 we would indeed have true singularities (that is, if we do not 
replace them by constant A and B solutions). In the former case as r — > rQ-\- we have 
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rU' ~ aU'^ , (116) 



where 



400 - 127V10 , 
a = ^ !^ 1.45 . (117) 



It then follows that 



40 - 13V10 



u ttV^' (lis) 

aln[r/ro) 

B ln(ln(r/ro)) , (119) 

a 



so that as r — > ro+ the curvature scalar diverges as 



2 



i?- (ln(r/ro))^^'"°^ , (120) 

and radial null geodesies terminate with finite affine parameter at r = tq. Similarly, in the 
latter case as r — > tq— we have 



rU' ~ (3U^ , (121) 



^ 400 + 127V10 , , 

/3= ^^ 9.88. (122) 



where 



40 + 13V10 
It then follows that 

ln(ln(ro/r)) , (124) 

so that as r ^ Tq— the curvature scalar diverges as 

(ln(ro/r))^('-^) . (125) 

and radial null geodesies are incomplete reaching r — rg within finite affine parameter. 

Thus, as we see, we have a sensible diluting non-singular solution, namely, the rj — —1 
solution with infinite volume and positive brane tension. In the next subsection we will 
discuss this solution in the context of the cosmological constant problem. 



D. The Brane Tension 

The last solution we discussed in the previous subsection, namely, that corresponding to 
the ^ > branch, is particularly interesting. This solution is diluting with positive brane 
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tension, so we would like to discuss it in a bit more detail. In fact, here we would like to 
give a general discussion of certain points relevant to the cosmological constant problem. 

Thus, we have been referring to A (or, more precisely, T = Mp~^A) as the brane tension. 
This quantity is indeed the tension of the brane E whose world- volume has the geometry 
of R'^"^'^ X Sf~^. Note, however, that a bulk observer at r > e does not see this brane 
tension - to such an observer the brane appears to be tensionless. Indeed, the warp factors 
are constant in the aforementioned solution at r > e. The non-vanishing (in fact, positive) 
tension of the brane S does not curve the space outside of the sphere Sf~^. Instead, it curves 
the space inside the sphere Sf"^, that is, at r < e. And this happens without producing 
any singularity at r < e, and with the non-compact part of the world-volume of the brane 
remaining fiat. 

Here it is important to note that the effective tension of the fat {d — l)-brane whose 
world- volume is R"'^^'^ (this {d— l)-brane is fat as it is extended in the extra {d— 1) angular 
dimensions) is also positive. It is not difficult to see that this brane tension is given by 

f = (A - (d - l)X)vd-iM^-^ , (126) 

where va-i — €^~^e'^'''~^^^'^^^Vd-i is the volume of the sphere Sf"-*^, va-i is the volume of the 
unit {d — l)-sphere, and we have taken into account that the radius of the sphere Sf~^ is 
not e but rather 

R = ee^(^) . (127) 

We therefore have: 

f^{d-2){K-{d- l))vd-iR'^-^M^~^ . (128) 

Thus, for K > 2 (which is part of the parameter space for the aforementioned solution) this 
effective fat brane tension is positive. 

Next, the d-dimensional Planck scale Mp is given by 

m|-2 ^ Vd-iM^~^ = Vd-iR'^-^M^-^ . (129) 

In particular, if we consider a 3-brane in 7D bulk, we have: 

A = 47r{K - 2)M| , (130) 
M| = 4rfM| , (131) 

where A is the 3-brane tension, Mp is the 4-dimensional Planck scale, Mp is the 6- 
dimensional Planck scale, and R is the radius of the extra 2-sphere (recall that we actually 
have a 5-brane with two dimensions curled up into a 2-sphere of radius R, while the radial 
direction transverse to this 5-brane is non-compact and has infinite- volume). 

A priori we can reproduce the 4-dimensional Planck scale Mp ~ 10^^ GeV by choosing 
R between R ~ millimeter and R ~ 1/Mp. The 6-dimensional Planck scale then ranges 
between^ Mp ~ TeV and Mp ~ Mp. On the other hand, it is reasonable to assume that 



^Note a similarity with [9]. 
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the analyses of [4] should give the same "see-saw" modification of gravity in the present 
case, in particular, once we take into account higher curvature terms on the brane. to obtain 
4-dimensional gravity up to the distance scales of order of the present Hubble size, we must 
assume that the "fundamental" 7-dimensional Planck scale^ 

Mp ~ (millimeter)-^ . (132) 

Let us see what range of values we can expect for the 5-brane tension T. 
Thus, the 5-brane tension 

f = AMj, = kRT'^MI ~ R-^Ml . (133) 

If we assume that the Standard Model fields come from a 6-dimensional 5-brane theory 
compactified on the 2-sphere^, then we might need to require that R <^ (TeV)~^. Then Mp 
ranges between 10^ — 10^ TeV and Mp, while the 5-brane tension ranges between (10^ TeV)^ 
and Mp. Note that a priori this is not in confiict with having the supersymmetry breaking 
scale in the TeV range. 

In principle, the above scenario a priori does not seem to be inconsistent modulo the 
fact that we still need to explain why the 6-dimensional Plank scale Mp is many orders 
of magnitude (30 in the extreme case where ~ Mp ~ Mp) higher than the seven- 
dimensional Planck scale Mp (see the last section for some speculations on this point). Note, 
however, that the same issue is present in any theory with infinite-volume extra dimensions. 
The diluting property of our non-singular solution could be suggestive as far as addressing 
the cosmological constant problem is concerned. In particular, we do have a solution where 
the non-compact 4-dimensional part of the brane is flat, yet its tension can take positive 
values in a continuous range. In particular, no flne tuning appears to be required in our 
solution. 

Before we end this section, let us see what kind of values of U (e— ) we would need to 
have in order to obtain a solution satisfying the above phenomenological considerations. 



^Actually, wc obtained the above solution by neglecting the Einstein-Hilbcrt term in the bulk. 
Here wc will simply assume that the aforementioned conclusion holds, and return to this issue 
when we study the more realistic model with both Einstein-Hilbert and Gauss-Bonnet terms in 
the bulk. 

■^Note that a priori we could attempt to have the Standard Model fields living on a 3-branc inside 
of the 5-brane (or at an orbifold fixed locus). This, however, would generically spoil the diluting 
property of the solution. Indeed, the aforementioned 3-brane inside of the 5-brane would have 
its own brane tension associated with it, which would generically have to be fine tuned. Another 
point worth mentioning here is that compactification on a 2-sphere of radius R <C (TeV)"^ need 
not in general result in supersymmetry breaking at a scale of order 1/R. Thus, one could imagine 
embedding our scenario in a higher dimensional theory (with perhaps some extra dimensions com- 
pactified) where the 2-sphere is fibered over a 4-dimensional base (such that the resulting space is 
supersymmetric;. in such 11-dimensional setup we might need, say, an additional Z2 orbifold to 
obtain a chiral 4-dimensional theory). Alternatively, the 5-brane could be a brane wrapped on a 
corresponding to an orbifold blow-up. 
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First, wc will assume that the Gauss-Bonnet parameter ^ ~ Mp^ (its "natural" value - see 
footnote 6). Then from the matching conditions we would obtain (C/(e— ) ^ 1): 

KRL/i ^ 9295 U^{e-) , (134) 

RL/^ ^ 4323 U^(e-) , (135) 

which gives «; ~ 2.15, and 

C/(e-) ~ lO''^ . (136) 

Here for definiteness we have assumed the extreme case ~ Mp ~ Mp, where L = 
Mp/M^ ~ lO^^'^ mm. This implies that e is very close to the would-be singularity Tq: 

^ - 1 ~ lO-^^t^ , (137) 

Note, however, that the warp factor 

B{e) ~ 7 , (138) 
which is due to the double logarithmic suppression in (124). 



IV. EINSTEIN-HILBERT AND GAUSS-BONNET TERMS IN THE BULK 

In this section we discuss a more realistic setup where we have both the Einstein-Hilbert 
as well as Gauss-Bonnet terms in the bulk. This substantially increases the complexity of 
the equations of motions we need to analyze, so we warm up with the example of a string 
in 5D bulk, and then turn to the most interesting case of a 3-brane in 7D bulk. 



A. String in 5D Bulk 

In this subsection we study the case of a string in five-dimensional bulk in the presence of 
both the Einstein-Hilbert and Gauss-Bonnet terms in the bulk action. The (rr), (rr) — (a/?) 
and {iiu) equations read, respectively: 

2 

{A' f + 4:N + S- 4^^ {2N^ + A'^S) = , (139) 



r 



2 



2{N-M) + S-P- 4C— 2N {N-M) + A'^ (S - P) = -e^AL5(r - e) , (140) 



1 



.B' 



2 



^2 1 ,^ s ^ 

M + 2N + 2P + S - 4^— {MS + 2NP) = --e^ (A - 2A) L5{r - e) . (141) 



A trivial solution for r 7^ e is given by A' = B' = 0. To study non-trivial solutions of this 
system, we define V = rA' and z = p^/^ (so that 2U = rz'/z, where we are using the 
definitions for U and p from the previous section), in terms of which (139) becomes: 

zV'^ + zU^ + AzVU - 12V^U^ + 4V^-z^0 . (142) 
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Thus, non-trivial solutions satisfy 

-2zU + CV3^2[/2 + i2zU^ _ i6^C/2 + ^2 + 4^ 
^ = + 4 ' 

AT . ^ , (144) 

M = (ry + - yC/) , (145) 

where ^ = ±1, and (140) reads 

2{z - 4UV){2UV -V^- rV) + {z - AV^){U^ - 1 - rU') = ]-ze^e^lL5{r - e) . (146) 

Away from the brane this gives a first order differential equation for [/ as a function of r, 
namely: 

rC/' = ^ , (147) 
or treating U as function of z we have: 

2.f/f^ = f , (148) 
dz /2 

where we have defined 

/i = mzU'^V - 2A(}zV^U^ + Az'^V^U + 34z'^U^V + 24zU'^V^ + 32V^ - 8zV + 
ASzUV^ + Sz'^U^ - Sz'^U - 2z^V - 256U^V^ + ASOV^U* - 8zV^ - Az'^V^ - 
192U^V^ + - AOzU^V + 16zUV^ , (149) 

/2 = GAzV^U - 6z^V - 96[/V3 - SzU'^V + 8zV - 32V^ - 8Vh , (150) 

for later convenience. 

Note that solutions must satisfy the following jump conditions implied by (146) and the 
(/iz/) equation of motion: 

[-z{2V + U) + 4V^u£_=^, (151) 



2C ' 

2e 



[4y(C/2 - 1) - z{2U + V)] ^{k-2)^ . (152) 



Let us rewrite these conditions as follows: 



[WU{V -U)+4V + z{U - V)] \t-^ ^ (153) 
%V'^U - 2z{U + 2V)] [_ = —K , (154) 

so that (153) only involves the coupling ^, whereas the ratio of (154) and (153) only depends 
upon K. We will consider two types of solutions: those that have constant warp factors A 
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and B for r < e, and non-trivial A and B for r > e - the exterior solutions; and those that 
have non-trivial A and B for r < e, and constant A and i? for r > e - interior solutions. 
Through the definition of the functions U and V we can rewrite (153) as follows: 

[4VU{V -U)+4V + z{U -V - 1)] = ±^ , (155) 

with the plus sign corresponding to exterior solutions and the minus sign corresponding to 
interior ones. Fig.5 shows the regions in the U — z plane where the function 4:VU{V — 
U) + 4V + z{U — V — 1) is positive (suitable for matching exterior solutions with ^ > 
or interior solutions with ,^ < 0) and negative (suitable for the other two possibilities) for 
the C = +1 case. We will now study the solutions with positive Gauss-Bonnet coupling ^ 
{i.e., z > 0) since this was the case that allowed non-singular solutions in the absence of the 
Einstein-Hilbert term in the bulk action. 



The ( = +1 Solutions 

In order to have exterior solutions we must have either U{e+) > 1 or U{e+) less than 
the negative zero of 4:VU{V — U) + AV -\- z{U — V — 1) shown in Fig.5. In the former 
case K is strictly negative, while k, is positive in the latter case. In both allowed regions 
— 1 < fi/{2zUf2) < 0. This means that the negative k solutions flow to lower values of U 
while z increases as r — > oo. A numerical study (using Maple 6) shows that as r increases, 
the solutions approach the region z^ 1,U >1. In this region (148) becomes (to the leading 
order in 1/z), 

1, (156) 



U -Idz 



which is integrated to give 



[/ ~ 1 + , (157) 



where C is a positive integration constant. This result is consistent with the r — > oo limit 
of 2Uz = rz', which can be integrated to give the asymptotic behavior of the warp factors: 
B{r) and A{r) approach constant values as r — > oo. Thus, these solutions are smooth and 
the volume is infinite. 

On the other hand, the positive k solutions flow towards higher values of U while z 
decreases as r — > oo. The numerical solution approaches the region z <^ 1, U < —1. By 
expanding 2zUf2/fi around 2; = we flnd the following asymptotic form of (148): 

16C/3 dU , , 

1 , (158) 



U'^-ldz 
which has the following solution, 

z 8([/2 + ln([/2 - 1)) + constant . (159) 
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Thus, a.s z 0+, U ^ < —1, a constant value. This in tiirn imphes that in this hmit 
(r — i> oo), B{r) ~ {U^ — l)ln(r/ro) and A{r) approaches a constant value. It is important 
to note that for these solutions the volume is finite. 

Let us now consider the possibility of interior solutions in which the volume is automati- 
cally infinite. For these to satisfy (155) with 2; > 0, U{e—) must he in the region in between 
solid lines in Fig.5 (the dash-dotted line separates the positive k (right) region from the 
negative k (left) one). These sohitions have potential finite distance singularities. Indeed, 
with initial conditions in the allowed region, the solutions flow to negative values of U and 
large values of z approaching the curve z — 2U^ as r decreases. A large U approximation of 
(147) along this curve yields, 

rU' - . (160) 

Thus, 

ln(r/ro) ' ^ ^ 

where Tq is an integration constant and U diverges as r ^ 1^0+, i-^-, we find a candidate 
for a finite distance singularity. However, since the asymptotic behavior of the warp factors 
implied by (161) is (up to overall numerical factors) 

A{r) ~ B{r) ln(ln(r/ro)) , (162) 

in the r — > ro+ limit, this singularity is just a coordinate singularity; it is indeed not difficult 
to check that in this limit radial geodesies are in fact complete. 

The ( = —1 Solutions 

Let us first study the exterior solutions. In order to have exterior solutions U (e+) must 
lie to the right of the solid line of Fig.6. In this region the string tension is negative. These 
solutions fiow towards the curve z = 2U^ for large values of U, and along this curve to the 
leading order in 1/U, 

rU' ~ . (163) 

Thus, 

ln(ro/r) ' ^ ^ 

where vq is an integration constant. U diverges in the limit r — > tq— and the asymptotic 
behavior of the warp factors in the same limit is as follows 

B ln(ln(ro/r)) . (165) 

However, radial geodesies are complete in the limit r tq—. And the volume of the extra- 
space is infinite. 
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Let us now consider the interior solutions. With 2;(e) > and U (e— ) to the left of the 
solid line of Fig. 6 we find solutions - with positive k, - that flow towards large values of z 
and U ~ — 1 as r decreases. In this region (148) takes the asymptotic form 

2z dU , , 

1 , (166) 



U + ldz 
which yields, 

-21n \U+l\r^ \n{z) + constant . (167) 

This implies that U ^ — 1 as 2; — +00. This is consistent with the r ^ limit; indeed, 
through 2zU = rz' we obtain, in this asymptotic domain, 

zr.1. (168) 

In this way we find the behavior of the warp factors in the r ^ limit: 

A ~ constant , (169) 
B ~ -21n(r/ro) . (170) 

These are infinite volume smooth solutions with positive k. 

Thus, as we see, in the case of a string in 5D bulk we have sensible smooth infinite- 
volume solutions with positive brane tension if in the bulk action we include both the 
Einstein-Hilbert and Gauss-Bonnet terms. Moreover, these solutions are diluting. 



B. 3-Brane in 7D Bulk 



In this subsection we study solutions with a fiat 3-brane in 7D bulk space, which is the 
most interesting case from the phenomenological point of view. The (rr), {jj.i') and {rr — a(3) 
equations of motion respectively read: 

6{A'f + 8N + S -12^^!^ {A'Y + + 8{A'fN + 2{A'yS^ = , (171) 



3{A'y + 3M + 6N + 2P + S - 



12^ ^{2^^ + {A'f{M + 2N + 2P + S) + AMN + + 2NP^ + 

-e^[A-2AjL5(r-e) = , 
2 

A{N -M) + S-P- 24^ ^|2Ar(Ar - M) + {A'f{2{N - M) + - P)| 



^e^AL(5(r - e) = . 



(172) 



(173) 



Using the notations of the previous subsection, the (rr) equation can be rewritten as 

C/2(72y2 _ _ 2UAV{z - 12V^) - z{6V^ - 1) - 12^^(2 -V^)^0 , (174) 



26 



which imphes 



^ ^ 72y^^^ ' 

where^ ?7 = ±1. 

Similarly, we can rewrite (173) as follows 

a(^z- 12V{U + V)^ {2U - V)V - rV + (z - 2W'^){U'^ - 1 - rU') = ^eRKL5{r - e) . 

(176) 

Thus, away from the brane we obtain the following first order differential equation for V: 



rV = 



9i 

92 



or treating V as a function of z we have: 



2zU 



dV _ 9i 
dz g2 ' 



(177) 



(178) 



where 



= 8Vh'^ - 192V^z - 4608V^U'^ + 576V^U - A320V^U^ - z^U^ + z^U + 2zV - 
72zV^ + mW^U + 1152\/^ + 576\^^ + Q72V^zU^ - IWz^U'^ - UV^z'^U + 
192V^zU + WSV^Uh - GOzV^U , 
c/2 = -12zVU^ + UAzV^U + lOS^y^ + 12zV - 5zV - 9>UV^U^ - 1152y^C/ 
-288^^ - 864^^ . 

To arrive this result we have replaced U' by the expression obtained by differentiating (174) 
w.r.t. r and we used the relation 



rz' = 2zU , 



(179) 



which follows from the definitions of U and z. 

The full solutions must satisfy the jump conditions coming from (176) and (172). They 
read, respectively. 



AiAV"^ - z)V - zU + 2AV^U _^ —RkL , 

3Vz + 2zU - AV^ - 2AV^U - 12V{U^ - 1)] = " '*)^^ 



(180) 
(181) 



In the limit ^ ^ oo, where we can neglect the bulk Einstein-Hilbert term, these solutions corre- 
spond, respectively, to the rj = ±1 3-brane solutions found in the previous section. 
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Let us rewrite these conditions as follows: 



12V^ + z{U -V)- 12VU^ 



32V' + 48V'U - 2zU - 8zV 



RL 
RL 

— —l—K 



(182) 
(183) 



As in the previous section we will consider interior {V{r > e) = 0, U{r > e) = 1) as well as 
exterior {V{r < e) — 0, U{r < e) — 1) solutions. For these, the above matching conditions 
take the following form: 



12V^ + z{U -V -1)- 12V{U^ - 1) 



32V^ + 48V^U - 2z{U - 1) - 8zV 



e± 



RL 



, RL 

±——K 



(184) 
(185) 



where the plus sign corresponds to the exterior solutions while the minus sign corresponds 
to the interior ones. 

In order to have non-trivial solutions consistent with the matching conditions, V{€+) 
(for exterior solutions), V{e—) (for interior solutions) and z{e) must be chosen in certain 
regions of the V — z plane which are identified in Fig. 7 and Fig. 8 for the 77 = +1 and rj = —1 
cases, respectively. In what follows we will study properties of the non-trivial parts of the 
solutions with initial conditions {V{e±), ^(e))in the allowed regions. 

The T) — +1 Solutions 

Let us first consider positive Gauss-Bonnet coupling, ^ > {i.e., z > 0). Exterior 
solutions are consistent with (184) if V{e+), z{e) lie in regions /, IV or V of Fig. 7 (through 
(185) regions / and V are consistent with positive k, while region IV is consistent with 
negative k). The solutions with initial conditions in region / flow to— V^^l,2;^l where 
(178) becomes 



dV V 

I— r 

az z 



(186) 



where a — 3(4 -|- vl0)/(5 -|- 14vl0)- In this approximation, 

z - C{-VY , U ~ -hV , 



(187) 

where h — {A -\- •\/T0)/6 and C is an integration constant. Now we can integrate (179) to 
obtain 



V 



1 



261n(ro/r) ' 



where ro > e is an integration constant and > —00 asr— >ro— (and z — > 
the asymptotic behavior of the warp factors in the r — > tq— limit is as follows: 



(188) 
+00). Thus, 
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^~^ln(ln(ro/r)), (189) 
C 

ex.p(B) r (ln(ro/r))"^ , (190) 

r 

where C is a constant. It is not difficult to show that radial null geodesies terminate with 
finite affine parameter as r — > tq— and the curvature scalar diverges in this limit implying 
that these solutions have a finite distance singularity^. 

The solutions with initial conditions in regions IV and V fiow to \/ ^ 1, 2; S> 1 ap- 
proaching the curve z = 12V'^ + 6/5 as r increases. In this approximation [/ ~ and (177) 
can be integrated to give, 

ln(ro/r) ' ^ ^ 

where ro > e is an integration constant and V diverges as r — > tq— . In this limit the warp 
factors have the following asymptotic behavior: 

A ~ -ln(ln(ro/r)) , (192) 

exp(S) ~ J (ln(ro/r)) , (193) 

where C is an integration constant. As in the previous case, these solutions have finite 
distance singularities. 

The interior solutions compatible with (184) have initial conditions in regions // and 
/// of Fig. 7 (the former is compatible with (185) for k, < while the latter is compatible 
with K, > 0). The solutions with initial conditions such that z{e) < 2AV^{e~) within region 
// fiow to the origin of the V — z plane. For z <^1, —V<^1 (178) becomes 

2^ . ^ . (194) 

dz z 

and in this approximation 



U ^U, = ^{24-a)/{72-a) , (195) 

where < a < 24 is an integration constant. The behavior of the warp factors as r ^ 0+ 
is: 

A ~ constant , (196) 
B^{U,- l)ln(r/ro) , (197) 

where ro is a constant. 

The solutions with initial conditions in the rest of region II flow to —V <^1, z ^ 1. In 
this approximation [/ ~ — 1 and the warp factors have the following asymptotic behavior: 



^This is the case for all the finite distance singularities that we mention throughout this section. 
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A ~ constant , (198) 
B ~ -21n(r/ro) , (199) 

where Tq is an integration constant. The r — singularity here is a coordinate singularity, 
in fact, it is not difficult to show that radial geodesies are complete in the r — > 0+ limit. 

Some of the solutions with initial conditions in region /// flow to the region where the 
solutions become complex (these solutions are imphysical), some flow to the origin of the 
V — z plane through region V, others flow to a flnite value of V = > as z ^ 0+, and 
there are solutions that flow to F ^ 1, 2; ^ 1. For the latter solutions in the large V, large 
z domain (178) becomes 

- - , (200) 

az z 

where a = 3(4 — vT0)vT0/(14vT0 — 5). In this approximation, U ~ —bV, where b = 
(4 — \/l0)/6, and we can integrate (179). We obtain: 

where Tq is an integration constant. V diverges as r ^ ro+ and the asymptotic behavior of 
the warp factors in this limit is given by: 

A~^ln(ln(r/ro)), (202) 

exp(S)~^(ln(ro/r))-t , (203) 

where C is a constant. These solutions have a flnite distance singularity. 

The solutions that flow to ~ 0.25 are also singular but in the r — > limit. In this 
limit for these solutions the warp factors are given by: 

A - -VM{ro/r) , (204) 
B - (1 - C/.)ln(r/ro) , (205) 

where C/* ~ 0.43 and Tq is an integration constant. 

The solutions that flow to the origin are also singular in the r — > limit. The analysis 
for these solutions is similar to that for the interior solutions with initial conditions in region 
// that flow to the origin. 

Let us now consider the ^ < case {i.e., z < 0). Interior solutions are consistent with 
(184) if (y(e— ), z{e)) lie in regions 7 or y of Fig. 7 which is consistent with (185) if k > 0. 
The latter solutions have the same properties as the interior solutions with initial conditions 
in region /// discussed above, while the former ones flow to the origin of the V — z plane 
and they are singular in the r limit. The analysis is similar to that for the interior 
solutions of region /// with the only difference that the constant a is negative in the present 
case. 

Exterior solutions are consistent with (184) if (y{e—), z{e)) lie in region VI of Fig. 7. If 
z{e) <^ 172.7\^^(e+) the solutions flow to —z :$> 1, —V -C 1. In this domain (178) becomes 
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2 dV V , , 

■=- . 206 

10 dz z ^ ^ 



In this approximation, [/ ~ 1, and thus A and B approach constant values in the large r 
Umit. Let us mention that these solutions are only consistent with negative values of k. 
For initial conditions in the rest of region VI the solutions have finite distance singularities. 
They flow to —V 3> 1, —z ^ 1 domain in which 

where a = 2(4 + VTO)/VW. We obtain 

where ro is an integration constant. V diverges as r — > tq— and the warp factors in this 
limit are given by: 

Ar.J=\n{\n(ro/r)) , (209) 

exp(B) - - (ln(ro/r))"^ . (210) 
r 

It is not difficult to show that here we indeed have a singularity at r = Tq. 

The T] — —1 Solutions 

We will first focus on the positive Gauss-Bonnet coupling case, ^ > {i.e., z > 0). 
Interior solutions are consistent with (184) if (y(e— ), 2;(e)) he in regions /, //, V or VI of 
Fig. 8 (/ and V are consistent with k > 0, while II and VI are consistent with k < 0). The 
solutions from regions / and // fiow to ^ ^ 1, —V ^ 1 approaching the curve z = 12V^. 
In this domain (177) becomes 

rV ~ V^ (211) 

with U ^V. We obtain 

ln(r/ro) ' ^ ^ 

where tq < e is an integration constant. V diverges as r — > rQ+ and in this limit the warp 
factors are given by: 

A~-ln(ln(r/ro)) , (213) 

where C is a constant. The above expressions are divergent in the r — > rQ+ limit, but the 
singularity is just a coordinate singularity. Both the affine parameter of radial null geodesies 
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and the proper time of time- like geodesies diverge in the same hmit as l/ln(r/ro). Thus, 
space is complete if we cut it at r = Tq. The values of k, consistent with (185) for initial 
conditions in region / range from along the dot-dashed curve growing to 2.15 for points 
away from it, while for initial conditions in region II k ranges from — oo (for z — 12V^, 
K ~ lOy/3 for —V » 1) to (along the dot-dashed curve). The solutions with positive k 
are of particular interest and we will come to them later. 

The solutions from regions V and VI flow to ^ ^ 1, \^ ^ 1 where (178) becomes 

~ ^ . (215) 

dz z 

where a — 3(4 -|- -\/IO)-\/IO/ (5 -|- 14-\/l0). In this approximation, 

z - , U - -hV , (216) 

where h — {A -\- -\/l0)/6 and C is an integration constant. Now we can integrate (179) to 
obtain 

— ^ — , (217) 

2Mn(r/ro) ' ^ ' 

where ro < e is an integration constant and V — > -|-oo as r — > tq-I-. The asymptotic behavior 
of the warp factors in the r — > ro-|- limit is as follows: 

A~^ln(ln(r/ro)) , (218) 

exp(S) ~ ^ (ln(r/ro))-§ , (219) 

where C is a constant. These solutions have a finite distance singularity. 

Exterior solutions are consistent with (184) if V{t+), z{e) lie in regions /// or IV . The 
former are consistent with k > while the latter are consistent with k < 0. The positive 
K solutions flow to the regions where the solutions become complex (these solutions are 
unphysical), while the negative k solutions flow to— y^l,2;^l. In this domain 

dV V , , 

-2- , 220 

dz z 

which gives V ~ C/z'^ for constant C. In this approximation ~ 1, which allows us to 
integrate (179) to obtain the behavior of the warp factors: both A and B approach constant 
values in the r —>■ oo limit. Let us mention that the volume for these smooth solutions 
is infinite. For these, the values of k range from approximately —0.77 to (along the 
boundaries of the region). 

Let us now study the negative Gauss-Bonnet coupling case {i.e., z < 0). V{e+), z{e) in 
regions 7 or y of Fig. 8 give exterior solutions consistent with (184) and (185) for positive 
values of k. Both have finite distance singularities. Those from region I flow to —V » 1, 
—z^l where (178) becomes 

~ ^ . (221) 

dz z 
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where a — 3(-\/l0 — 4)-\/l0/ (5 — 4-\/l0). In this approximation, 



z - C{-VY , U - -bV , (222) 
where b — {A — ^/W)/6 and C is an integration constant. We obtain 

V ~ - — , (223) 

2Mn(ro/r) ' ^ ^ 

where ro > e is an integration constant and V— > — ooasr— >ro— . The asymptotic behavior 
of the warp factors in the r — > tq— hmit is as follows: 

A ^ln(ln(ro/r)) , (224) 

C 

exp(B) (ln(ro/r))"^ , (225) 

r 

where (7 is a constant. 

The solutions from region V flow to the origin of the V — z plane. For V <^ 1, —z<^ 1, 

2^ ~ . (226) 

dz z 



Thus, z ~ — aV^, where a > is an integration constant, and U — (24 + a)/ (72 + a) . 
Integrating (179) we obtain 

V ~ -Cr-^* , (227) 

where C is an integration constant. V diverges as r ^ 0+. The warp factors are also 
divergent in this limit. 

Interior solutions are consistent with (184) if V{e—), z{e) lie in region VIII of Fig. 8. 
Throughout this region only negative values of k are consistent with (185). For V{e—) < 
the solutions flow to —V <^ 1, —z ^ 1 and in this domain we have U ~ —1. The behavior 
of the warp factors in the r ^ limit is given by: 

A ~ constant , (228) 
B ~ -21n(r/ro) . (229) 

The singularity in this case is just a coordinate singularity as geodesies are complete in the 
r — > 0+ limit. 

For V{e-) > with z{e) ^ - l?>AlV^{e-) the solutions flow to F < 1, > 1 with the 
same asymptotic behavior of the warp factors as in the previous case. For initial conditions 
in the rest of region VIII the solutions flow to y ^ 1, —z <^ 1 with the same behavior 
as for the interior solutions of region V for the positive z case. They have finite distance 
singularities. 
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Diluting Solutions with Positive Tension 



Let us now come back to the ?] = — 1, { > interior solutions with initial conditions 
in region I of Fig.8. Prom our discussion in the previous section, these solutions are of 
particular interest as some of these solutions have positive 3-brane tension (the ones with 

K > 2) and, furthermore, they are diluting. 

We can for instance take the ^ <^ |y| ^ 1) case in which (184) and (185) become, 
respectively, 

RL/^ ^ -11.77^3(6-) , (230) 
kRL/^ ^ -25.3V\e-) , (231) 

which gives k ^ 2.15, and |y(e— )| fsi 10^°, and e is very close to the would-be singularity r* 
(here for definiteness we have assumed the extreme case ~ Mp ~ Mp) : 

^ - 1 - 10-3° . (232) 

Note that the singularity at r = r* > e would be there if we took the interior solution 
for r < e and continued it for values r > e. However, in our solution (just as in the case 
without the Einstein- Hilbert term) there is no singularity as for r > e the warp factors are 
actually constant, and this is consistent with the matching conditions at r = e. In particular, 
the solutions are smooth everywhere, just as their counterparts from the previous section. 
Thus, as we see, in the presence of both the Einstein- Hilbert and Gauss-Bonnet terms in the 
bulk action we have sensible smooth solutions with positive brane tension. Moreover, these 
solutions are diluting, that is, they exist for a range of values of the brane tension (note that 
the Gauss-Bonnet coupling in these solutions is positive). 

Before we end this subsection, let us emphasize some important points. In the diluting 
solutions we just discussed, we have a coordinate (but not a real) singularity at some finite 
r — ro < e. As we mentioned above, this coordinate singularity is harmless as the corre- 
sponding geodesies are complete. Note that in the case without the Einstein-Hilbert term 
the corresponding coordinate singularity is at r = 0. The reason why is the following. If we 
start with a solution corresponding to both the Einstein-Hilbert and Gauss-Bonnet terms 
present in the bulk, to arrive at the solution corresponding to only the Gauss-Bonnet term 
present in the bulk we must take the limit ^ — > oo, Mp~^^ — fixed. It is then not difficult 
to check that in this limit the coordinate singularity at r = tq continuously moves to r = 0. 
Also note that, since the singularity at r = ro in the full solution is a coordinate singularity, 
we can consistently cut the space at r = tq. The geometry of the resulting solution then 
is as follows. In the extra three dimensions we have a radially symmetric solution where 
a 2-sphere is fibered over a semi- infinite line [ro,oo). The space is curved for Tq < r < e, 
at r = e we have a jump in the radial derivatives of the warp factors (because r = e is 
where the brane tension is localized), and for r > e the space is actually flat. So an outside 
observer located at r > e thinks that the brane is tcnsionlcss, while an observer inside of 
the sphere, that is, at r < e sees the space highly curved (and it would take this observer 
infinite time to reach the coordinate singularity at r = Tq). This is an important point, in 
particular, note that we did not find smooth exterior solutions where the space would be 
curved outside but flat inside. And this is just as well. Indeed, as was shown in [7] , if we have 
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only the Einstein-Hilbert term in the bulk, then we have no smooth solutions whatsoever 
(that is, smoothing out the 3-brane by making it into a 5-brane with two dimensions curled 
up into a 2-sphere does not suffice) . What is different in our solutions is that we have added 
higher curvature terms in the bulk, which we expected to smooth out some singularities. 
But serving as an ultra-violet cut-off higher curvature terms could only possibly smooth out 
a real singularity at r < e, not at r > e. And this is precisely what happens in our solution 
- the presence of higher curvature terms ensures that we have only a coordinate singularity 
at r < e instead of a real naked singularity as would be the case had we included only the 
Einstein-Hilbert term. 

C. Implications for the Cosmological Constant Problem 

In this section we saw that in the presence of the Einstein-Hilbert and Gauss-Bonnet 
terms in the bulk action we have smooth infinite-volume solutions which exist for a range 
of positive values of the brane tension (the diluting property). These solutions, therefore, 
provide examples of brane world scenarios where the brane world- volume can be flat without 
any fine-tuning or presence of singularities. Is this then a solution to the cosmological 
constant problem? 

To answer this question, we need to address some issues here. First, note that we have 
chosen a particular combination of quadratic in curvature terms in the bulk action, namely, 
the Gauss-Bonnet combination. One could argue that this is a fine-tuning as we have to fix 
two independent parameters at the quadratic level to obtain the Gauss-Bonnet combination 
(note that one out of three a priori independent parameters corresponds to the Gauss- 
Bonnet coupling, which is arbitrary in our solutions). However, we suspect (albeit we do 
not have a proof of this statement) that even for generic higher curvature terms solutions 
with the aforementioned properties should still exist. In particular, we suspect that the fact 
that we found non-singular solutions has to do with including higher curvature terms in the 
bulk rather than with their particular (Gauss- Bonnet) combination, which we have chosen 
to make computations tractable. 

A more serious issue here has a phenomenological origin. Thus, as we discussed in 
subsection D of section HI, the 6-dimensional Planck scale on the 5-brane (whose world- 
volume is a product of the 4-dimensional Minkowski space-time and a 2-sphere of a radius 
R) Mp » Mp - as was explained in [4], we must have the 7-dimensional bulk Planck 
scale Mp ~ (millimeter)"^, so that the four-dimensional laws of gravity persist up the 
the distance scales of order of the observable Hubble size. Naturally, here we can ask the 
following question: Why is the 6-dimensional Planck scale on the 5-brane many (up to 30 
in the extreme case where Mp ~ Mp) orders of magnitude larger than the 7-dimensional 
Planck scale? We would like to give one speculative scenario for obtaining such a hierarchy 
of Planck scales. Thus, let the 5-brane theory be a (non-conformal) large N gauge theory. 
Then the Planck scale induced on the brane due to quantum effects is expected to be of 
order of Mp ~ N'^Mp. The required rank of the gauge theory in this case would have to be 
up to A?" ~ 10^° (in the extreme case). It would be interesting to understand if one could 
accommodate the Standard Model in such a scenario. 

There are many interesting open questions to be addressed in scenarios with infinite- 
volume extra dimensions. As was originally pointed out in [6,10,1], these scenarios offer a new 
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arena for addressing the cosmological constant problem. Wc hope our results we presented 
in this paper are at least encouraging in this context. And addressing the aforementioned 
open questions definitely seems to be worthwhile. In this light, we would like to end this 
paper with the following poem by a 19th century Russian poet Yuri Lermontov (translation 
from Russian by Z. Kakushadze): 

The Sail 



A lonely sail seeming white, 
In misty haze mid blue sea. 
Be foreign gale seeking might? 
Why home bays did it flee? 

The sail's bending mast is creaking. 
The wind and waves blast ahead. 
It isn't happiness it's seeking. 
Nor is it happiness it's fled! 

Beneath are running azure streams. 
Above are shining golden beams. 
But wishing storms the sail seems. 
As if in storms is peace it deems. 
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FIGURES 
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FIG. 1. Solutions for A{r) (in e = 1 units with ^(e) set to zero) for various values of a = kL, 
b = (k — 2)L (and ^ - see below). From top to bottom: ^ = 1, a = 2, 6 = —1; ^ = 1, a = 1, 6 = — 1; 
^ = 1, a = 2, 6 = —1.2; ^ = 1, a = 1, 5 = —1.2 (The ^ = 2 counterparts of the lines for B{r) in 
Fig. 2 coincide with the ^ = 1 lines.) This plot shows that Solution A for the membrane is indeed 
diluting. 
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FIG. 2. Solutions for B{r) (in e = 1 units) for various values of a = kL, b = (k — 2)L and ^. 
From top to bottom: ^ = 2, a = 1, 6 = -1; ^ = 2, a = 1, 6 = -1.2; ^ = 1, a = 2, 6 = -1; ^ = 1, 
a = 1, 6 = -1; ^ = 1, a = 2, 6 = -1.2; ^ = 1, a = 1, 6 = -1.2. This plot shows that Solution A for 
the membrane is indeed diluting. 
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FIG. 3. Solutions for A(r) (in e = 1 units with ^(e) set to zero) for A = 1 (see below) and 
various values of a = L/^. Prom top to bottom: a = 10^, a = 10^, a = 10^ and a = 1. (The A 7^ 1 
counterparts of the solid lines for B{r) in Fig.4 cannot be distinguished from the a = A = 1 line.) 
This plot shows that the ^ > branch of the 77 = — 1 3-brane solution is indeed diluting. 
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FIG. 4. Solutions for B[r) (in e = 1 units) for various values of A and a = L/^. From top to 
bottom: o = 1 with A = .1, A = 1, A = 10, A = 100 (solid lines) ; and A = 1 with a = 10~^, 
a = 10^, a = 10^ (dashed lines) ; This plot shows that the ^ > branch of the rj = —1 3-brane 
solution is indeed diluting. 
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FIG. 5. The solid lines are zeros of AVU{V + U) + W + z{U - V - 1) while the dot-dashed 
lines are the zeros of —2z{2V + [/ — !) + SV'^U. In the region in between dashed lines there are no 
real solutions of (142). 
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FIG. 6. The solid line represents the zeros of AVU{V + U) + 4V + z{U - V - 1); inside the 
dashed pockets there are no real solutions of (142). 
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FIG. 7. The solid lines are zeros of 12V^ + z{U - V - I) - 12V{U'^ - 1) that approach the 
following curves for \z\ » 1: I — II boundary: V = —2/3, /// — IV boundary: z = 12V^ and 
I -VI boundary: z = -27.3^^. The dot-dashed line between II and /// is 2; = 72^^ while the 
one between IV and V that approaches z = 9.61^^ for 2; S> 1 consists of zeros of k. Inside the 
dashed pockets there are no real solutions of (174). 
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FIG. 8. The solid lines are the zeros of 12V^ + z{U -V -I)- 12V{U'^ - 1) that approach the 
followina; curves for \z\ > 1: II - III boundary: z = 12.7V^, V - VIII boundary: z = -3.66^^ 
and / — VIII boundary: V=-2/3. The dot-dashed line between /// and IV is z = 72V^ while the 
one between I and II representing the zeros of k approaches z = — 3.66^^ for V ^ 1 (the third 
dot-dashed line consists also of zeros of k). Inside the dashed pockets there are no real solutions 
of (174). 
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